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Tomography of many-body weak values: Mach-Zehnder interferometry 
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We propose and study a weak value (WV) protocol in the context of a solid state setup, specifically, 
an electronic Mach-Zehnder interferometer. This is the first specific proposal to measure both the 
real and imaginary part (i.e., complete tomography) of a WV. We also analyze the manifestation 
of many-body physics in the WV to be measured, including finite temperature and shot-noise-like 
contributions. 
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The problem of non-invasive measurement of a quan- 
tum system is interesting both from the view point of 
foundations of quantum mechanics and possible appli- 
cations (e.g. quantum computation). As opposed to the 
standard strong measurement procedure described by the 
projection postulate [l[ , weak measurement of an observ- 
able, while weakly disturbing the system, provides only 
partial information on the state of the latter. Normally 
the outcome of a weak measurement of an observable is 
its expectation value, which is the weighted average of 
the eigenvalues of this observable. Major step in formu- 
lating alternatives to standard measurement procedure 
was the proposal of a weak value protocol @ , consisting 
of a weak measurement (of A), followed by a strong one 
(of B). The outcome of the first is conditional on the 
result of the second (post-selection). Specifically, if the 
system is prepared (preselected) in the state \xi), and 
is postselectcd in the state |x/)> the weak value of the 
operator A is 



(xilMxi) _ (xi\n Xf A \xi) 



(Xf\Xi) 



(Xil^xAXi 



(1) 



where IT X/ = |x/} (xf\ is a projection operator. 

In general, the weak value is a complex number and 
(as opposed to a projective expectation value) its real 
part can be out of the range of the eigenvalues of A. In 
order to obtain non-trivial weak values (e.g., outside the 
range of the eigenvalues of A; complex; negative when A 
is positive definite), the preselected state should not be 
an eigenstate of the measured operator, and [A, B] ^ 0. 

WV may allow us to explore some fundamental aspects 
of quantum measurement, including access to simultane- 
ous measurement of non-commuting variables 

HH; de- 
phasing andphase recovery Q ; correlation between mea- 
surements [J; and even new horizons in metrology 0, 0] • 
While some aspects of WVs have been demonstrated in 
optics based setups 0] , the implementation of weak val- 
ues in the context of solid state physics is very new P.[lfj|. 

Weak values are not just an artifact of a convoluted 
definition, but indeed do emerge as the outcome of weak 



measurement of a pre- and post-selected states. Micro- 
scopically we consider the coupling of a system to a de- 
tector with a Hamiltonian H = H$ + H-q + H int . The 
interaction Hamiltonian is H mi = Xg(t)pA. Here p is 
the momentum canonically conjugate to the position of 
the detector's pointer, q, and Xn(t) (A < 1) is a time 
dependent coupling constant Following the weak 

measurement and the post-selection steps, the expecta- 
tion value of the coordinate of the pointer (initially equal 
to qo) is given by (q) = q — XRe[^,(A) ( f > ], i.e., the shift 
in the pointer is proportional to the real part of the WV. 
Under more general conditions (i.e., p, q are not canon- 
ically conjugated]the imaginary part of the WV may be 
meaningful too [l2j. An experimental procedure which 
will provide for a full "tomography" of (both the real and 
imaginary parts of) the weak value remains a challenge. 

Here we report the first systematic study of complex 
WVs in the context of many electron solid state system, 
proposing an experimental procedure which will provide 
for the full "tomography" of (both the real and imag- 
inary part of) WVs. Employing building blocks which 
are accessible within current technology [l3j . our pro- 
posed protocol and the ensuing predictions are amenable 
to experimental test. In particular, addressing a "sys- 
tem" and a "detector" which are represented by an elec- 
tronic Mach-Zehnder interferometer (MZI) [l3| (Fig. Q}, 
(i) We propose how to retrieve both the real and the 
imaginary parts; (ii) We show that the introduction of 
both a non-pure state and finite temperature (T) mod- 
ify the WV and reduce the visibility of Aharonov-Bohm 
(AB) oscillations; (iii) We show how many-body effects 
lead to reduction of the weak value's visibility, as a func- 
tion of voltage bias. In particular Im[ X/ (A) y . 1 main- 



tains its single particle value while Re[ X/ {A) x ] is mod- 
ified by an excess (non-equilibrium) noise and thermal 
noise terms; (iv) We show that for the system at hand 
the WV is related to a system-detector current-current 
cross correlator [T3 |. 

The chiral leads of the MZI are realized by the edge 
states of an integer quantum Hall setup. Throughout 
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FIG. 1: (a) Electronic MZI represented by edge states (full 
lines) in a Hall bar. Inter-edge tunneling (dashed lines) takes 
place at the quantum point contacts, (b) Scheme of a double 
electronic MZI. Electrons on edge states 2 and 3 interact via 
Coulomb interaction. 



this analysis we will consider non-interacting single edge 
channels {v = 1) [l^]. We do account here for the system- 
detector interaction and for the inherent many-fcrmion 
physics. 

Our model Hamiltonian for the system (S) and the 
detector (D) is [HI 



H s = -iv ^2 dx : iljl l (x)d x ip m (x) : 
m=l 

+I>1(0)V> 2 (0) + T h <ip\{h)il> 2 (h) + H.c. 
H D = -iv y^ldx : ip} n (x)d x ^ m (x) : 

m=3 J 

+r c v4(o)v> 4 (o) + r^UWMh) + h.c. 



(2) 



(3) 



Here ?/4(a 



1 v^ c 
Jr. Z—'k 



e ikx dl 



^-fc-oo c "fcm is thc creation op- 
erator at point x on arm m (the system's length, L, is 
assumed to be much larger than any other length scale 
in the problem). The T's are inter-edge tunneling am- 
plitudes, normal ordering is with respect to the vacuum 
(states are occupied for k < and empty otherwise), v is 
the edge velocity and h = 1. We assume a weak capaci- 
tive interaction between the system and the detector: 



where the support of U(x,y) is < x,y < min m (l m ). 
Without loss of generality we assume that sources 5*1 and 
53 are voltage biased (V and V3), and use the current at 
-D4 as the weak measurement's pointer. 

Hs is readily diagonalizcd. In thc Hciscnbcrg repre- 



sentation the system's ^ operators are 



ik{x—vt) 



bkrn 1 



m • 



x < 0; 

< X < l ra ] 
I m *£■ 



(4) 



tering matrices Si = 



(i = a, b) and, for 



where m = 1,2 and = S A a.k, = S B b^ with scat 

n t, 

-n r, 

equal arm lengths, n = \{2v) 2 - |r. i | 2 ]/[(2i;) 2 + |r,| 2 ], 
U = AivTi/[{2v) 2 + \T t \ 2 ]. The effect of the system's AB 
phase and unequal arm's lengths can be absorbed into 

' 2 . Here 6 includes the con- 



T A = T A e^ 2 , T B =T B e~* 
tributions of both the AB flux, and the orbital phase: 
4> = 4> + k(h - h), <\> = 2tt^ (T a ,T b > 0). A similar 
expression is available for Hr> as well. 

As a first step we consider the injection of a single 
particle (in the scattering state k) into the system 17 1 
(e.g., at the source SI). The preselected state is then 

4lo) 



S\) = a' kl |0). The post-selected state is I-D2) = 
(detection at D2). By definition (JTJ) the WV of the charge 

/e = 

(pure 



at the intermediate segment of arm i is D2{Qi)%\/ e 



-, yielding D2(Q 2 ) s \h 



t*A r B 



t\r B +e i <l>T A fi 



state). This quantity may assume complex values de- 
pending on the choice of parameters (transmission and 
reflection amplitudes and $). In particular, its real part 
may be larger than 1 or even negative, clearly outside 
the allowed range of the number operator expectation 
value. But, as expected, D2(Qi)|i +D2 (<32)|^ = e. The 
individual weak values are <E>o periodic with the flux. 

In a more general case, when the system is initially not 
in a pure state (e.g., finite temperature, entanglement, 
etc.), eq. ([I]) is not applicable, and we have to rewrite 
it in terms of the density matrix, p, of the initial state, 
yielding [§]: 



Xf 



Tr(Tl Xf Ap) = (xf\Ap\ X f) 
Tr ( u x f P) (Xf\p\Xf) 



D2{Q2lSl 



For example assume that a voltage bias V is applied to 
the source SI, which is held at temperature T. Then 

p = l/LY, k [f{vk - eV) - /(i>fc)]aj-i a >si> witn f( x ) the 
Fermi-Dirac distribution at zero voltage bias, and 

t* A r B {t A r B + K(V, T)e-^r A t B ) 
\r A t B \ 2 + \t A r B \ 2 + 2K(V,T)Re[e^r A t* B t A r B } 

(5) 

where v = eVM/{2v)+^ K(V, T) ^ ™ B ™£ B %$$> 
and Al = I2 — h- We observe that the broadening of the 
energy and non-pure nature of the initial state lead to 
the suppression of the interference term, which in turn 
leads to the suppression of the imaginary part of the WV 
and reduces the real part to a "standard" value lying in 
the interval [0, 1]. 

We next study the generalization of this protocol to 
a many-body (albeit non-interacting) system. The prese- 
lected state is given by the stationary voltage biased (V) 
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non-equilibrium state of the system. Post-selection will 
be defined by the detection of one electron at the drain 
D2, i.e. within [x — e, x + e], e — ► around x > 
The projection operator for this state is : i/jl(x)ip 2 (x) ■■ 
For a system with linear dispersion relation, the current 
operator is l2(x,t) = ev : ipl(x,t)ip 2 (x,t) :; it follows 
that we can employ I 2 /(ev) as a projection operator for 
the post-selection. We consider the WV of charge (or cur- 
rent) at the intermediate segment of arm 2. It is given by 



the operator Q 2 (xq) = ev fl^^dt ■ ipl(xo,t)ip 2 (xo,t) :, 

where r is an infrared cut-off, r ^> and < xq < l-x- 
Such an operator can be measured e.g. by adiabatically 
switching on (off) the coupling to the measurement de- 
vice at t < (t > 0) [18j |, From (fT]), and employing 
(aLaki) = f(vk - fii) with fit = eV, [i 2 = 0, we ob- 
tain ngi: 



in \ MB 

D2 \W2)si 



(h(x,0)Q 2 (x )) 
(h(x,0)) 



VLZ k (fev-f )A 



-[E zE^A-^(i-/ w )]+^E^ 



r 



f eV -fo)Bu ■ (6) 



Here Ay = (S' A S' B ) i2 (S B S 'a)2 j , % = (<5^) i2 (5A)2j and 
fij, stands for the Fermi distribution function f(vk — /i). 
The summation over k is replaced by energy integrals, 
l/LYjk ~* l/(27r«) JdE. At zero temperature the WV 
© is: 

D2{Q2)T=D2{Q2)f 1 +(N-l){Q 2 ) (T = 0). (7) 

Here N = is the number of electrons injected into 
the MZI during the measurement time. The first term on 
the r.h.s is the single particle WV §5§. The second term 
is proportional to the single particle expectation value of 



Q2 = \tA\ 2 - This can be interpreted in the following way: 
the post-selection enforces a detection of one electron at 
point x at time t. The correlation between Q 2 and a post- 
selection yields the first term. All other electrons which 
arc injected into the system during the measurement are 
not constrained by the postsclcction and therefore give 
rise to an additional contribution to U2 (0.2)^3^ ■ 



Consider now a finite 
write the weak value as 



temperature. 

K B ± J 2-kv dE 



Using 
we can 



J 



in \ MB 

D2 [W2/SI 



Re 



F(V,T) (\r B r A t A \ 2 - \t B t A r A \ 2 ) + (\r A \ 2 - \t A \ 2 ) K{V 7 T)e l <?r A t A r B t\ 



\r A t B \ 



-iJm 



D2 (Q2) S li 



\t A r B \ 2 + 2K(V,T)Re(e i n A r B r A t* B ) 
k B T\r B \ 2 



eV(\r A t B \ 



\t A r B \ 2 
I 



2K(V,T)Re(ein A r B r A t* B )) 



+ (Q 2 )(N-1). (8) 



Here F(V,T) = coth[eV/(2fc B T)] - 2k B T/(eVV 
K{V,T) is given by Eq. ©, and (2fl| 
K(V,T) = 2TTk B T[sm{eVAl/2v)coth{eV/2k B T) - 
k B TAl/vcos(eVAl/2v)]/[eVsmh{irk B TAl/v)}. This is 
the central result of the letter. The many body weak 
value in electronic MZI consists of four terms: (i) A 
term which is related to the real part of the single 
particle weak value corrected by finite temperature 
Fermi sea effects; (ii) The imaginary part of the single 
particle weak value; (iii) A contribution due to thermal 
fluctuations [2l[; (iv) Many body "excess noise". The 
dependence of the many-body WV on voltage and 
temperature is presented in Fig. [2] Note the asymmetry 
between the imaginary and the real parts of the WV: the 



former coincides with the imaginary part of the single 
particle weak independently on the voltage bias and 
temperature. When the temperature is smaller than the 
voltage bias, k B T <C eV the many-body VW reduces to 
the expression in Eq.([7]), while at high temperature the 
third term dominates and thermal fluctuations wash out 
the SP weak values. 

We now describe the measurement procedure lead- 
ing to the WVs. We should calculate the expectation 
value of the pointer's variable, which, in our case, is 
the time average of the outgoing current in the drain 
D4. The latter is given by the expectation value of 
i 4 (x') — l/r j_^, 2 dtli{x',i) (x' > U), in the detector's 
state after postsclcction. Post-selection is realized by act- 
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FIG. 2: Real and Imaginary part of the many-body WV 
as a function of the Aharanov-Bohm phase. Panels (a), 
(b) correspond to different plots for different values of tem- 
perature: TA/l/(hv) =0 (dotted), 0.05 (dash-dotted), 0.15 
(dashed), 0.5 (long-dashed), 1.2 (full). V = 5fieV. Pan- 
els (c), (d) correspond to different plots for different volt- 
age biases: VAl/(hv) =0.2 (dotted), 0.3 (dash-dotted), 0.5 
(dashed), 1 (long-dashed), 2 (full). T = 3mK. In all plots 
fiv/Al = 5 x 10~ 6 eV, r = 2 x 10" 10 s, t A = r A = 1/V2, 
t B = y/2r B = V2/VS. The range for which Re[_D2(Q2)gi] is 
beyond the interval of standard (strong) values is shadowed. 



ing on the system's state with the projection operator 
I 2 (x,t). Therefore the effective action of the detector 
after performing the post-selection and tracing over the 
system's state is 



JS 



(eff) 



(9) 



and 



7 (cff) 



function 
We 



is 
cal- 



the effective partition 

culate the expectation value of Ii(x',t) using 
this effective partition function, which yields 
h{x')) = l/4 eff) JV[^ 3 }V[^ i }h(x')e iS n t) = 

I4,(x')l2(x,t)\ I ll2(x,t)\. Hence, in order to obtain 

a weak value we consider the correlator (I 2 (t)l4(t')). 
Since the system is out of equilibrium we employ the 
Keldysh technique to calculate this correlator. After 
some algebra we obtain to first order in the interaction 



1 



t+r/2 



dt' {h{t)h{t')) = (h) (h)o 



H-t/2 

I™{( D2 (Q 2 )f*)( Di (Q 2 )f 3 B )}) , (10) 



2A 



where A 



1 J dxdyU(x,y), and 



refers to the A = 



expectation value. Indeed , dividing by (I 2 ), we see that 
the shift in the detector's current is proportional to the 
WV ((6]). Since here the "system" and the "detector" play 
a symmetric role, the result can be stated in a system- 
detector symmetric fashion: for weak enough interaction, 



the current-current correlator is proportional to the prod- 
uct of WV of the respective system and detector op- 
erators which appear in i?i n t. Notice that varying the 
Aharonov-Bohm flux in the " detector" MZI one can ex- 
plore the real and the imaginary parts of the "system's" 
WV. 

The analysis presented here is the first significant step 
towards a complete characterization of weak values in in- 
teracting systems. The visibility of flux sensitive many- 
body WVs is reduced by both T and V: The single 
particle WV is supplemented by thermal noise and by 
a voltage dependent many fermion term. This WV to- 
mography and the tuning of the respective noise terms 
are all amenable to experimental verification. 

We acknowledge illuminating discussions with Y. 
Aharonov, A. Stern, and L. Vaidmann. This work was 
supported by the Minerva Foundation, US-Israel BSF 
and the DFG project SPP 1285, and DFG Priority Pro- 
gramme " Semiconductor Spintronics" . 



[9 
[10 

[11 
[12 

[13 



Deceased 

J. von Neuman, Mathematische Grusndlagen der Quan- 
temachanik (Springier- Verlag, Berlin 1932) 
Y. Aharonov, D. Z. Albert, and L. Vaidman, Phys. Rev. 
Lett. 60, 1351 (1988); Y. Aharonov and L. Vaidman, 
Phys. Rev. A, 41, 11 (1990). 

W. Hongduo and Y. V. Nazarov, |cond-mat/0703344 
(2007). 

A. Jordan and M. Buttiker, Phys. Rev. Lett. 95, 
220401(2005). 

I. Neder, F. Marquardt, M. Heiblum, D. Mahalu, and 
V.Umansky, Nature Phys. 3, 534 (2007); I. Neder, 
M. Heiblum, D. Mahalu, and V. Umansky, Phys. Rev. 
Lett. 98, 036803 (2007). 

E. Sukhorukov, A. Jordan, S. Gustavsson, R. Leturcq, 
T. Ihn and K.Ensslin, Nature Phys. 3, 243 (2007); A. 
Di Lorenzo and Y. V. Nazarov, Phys. Rev. Lett. 93, 
046601(2004). 

O. Hosten, and P. Kwiat, Science 319, 787 (2008) 

G. J. Pryde, J. L. O'Brien, A. G. White, T. C. Ralph, 

and H. M. Wiseman, Phys. Rev. Lett. 94, 220405 (2005); 

A. M. Steinberg, Phys. Rev. Lett. 74, 2405 (1995); H. M. 

Wiseman, Phys. Rev. A 65, 032111 (2002). 

Alessandro Romito, Yuval Gefen, and Yaroslav 

M. Blanter, Phys. Rev. Lett. 100, 056801 (2008). 

N. S. Williams and A. N. Jordan, Phys. Rev. Lett. 100, 

026804 (2008). 

In this example we have assumed that the free Hamil- 
tonians of the system and the detector vanish and that 
g(t)=6(t-t ). 

See e.g., A. M. Steinberg, Phys. Rev. A 52, 32 (1995); L. 
M. Johansen, Phys. Lett. A 322, 298 (2004); R. Jozsa, 
Phys. Rev. A 76, 044103 (2007); We have learned of 
the interesting work of A. Di Lorenzo and J. C. Egues, 
|arXiv:080lT814] . 

Y. Ji, Y. Chung, D. Sprinzak, M. Heiblum, D. Mahalu 
and H. Shtrikman, Nature 422, 415 (2003); P. Roul- 
leau, F. Portier, D. C. Glattli, P. Roche, A. Cavanna, G. 



5 



Faini, U. Gennser, and D. Mailly, arXiv:0704.0746v2]; I. 
Neder, N. Ofek, Y. Chung, M. Heiblum, D. Mahalu, and 
V. Umansky, Nature 448, 333 (2007). 

[14] We believe that a qualitatively similar statement (con- 
cerning a system-detector two operator cross-correlator) 
applies to weak value protocols at large. 

[15] While this simplified picture may not account for some 
aspects of the problem [dephasing (cf. F. Marquardt and 
C. Bruder, Phys. Rev. Lett. 92, 056805 (2004)), "lobe" 
structure, etc.], we find it as a relevant reference model 
to account for the weak value physics 

[16] J. T. Chalker, Yuval Gefen, and M. Y. Veillette, Phys. 
Rev. B 76, 085320 (2007) 

[17] Note that a single particle MZI can be mapped onto the 
problem of spin 1/2 particle. Indeed, choose the z-axis 
such that |T). = b[\0) and |J.) = b\ |0). Then a z = 



b\bi — &2&2, the scattering matrices rotate the spin in, 
say, the xz plane and the AB flux rotate the spin around 
the z-axis. 

[18] r 2> L/v allows us to replace the weak measurement 
signal taken over a finite segment by an operator defined 
at a point (xo). 

[19] This expression is readily generalized to the situation in 
which bias is applied to the source S a , the post-selected 
electron is detected at Dp, and the weakly measured op- 
erator is Q 7 . In this case one has to replace the index 2 
in the definition of Aij by f), and in Bij by 7. The index 
1 in /1, An, and Bn is replaced by the index a. 

[20] V. S. -W. Chung, P. Samuelsson, M. Buttiker, Phys. Rev. 
B 72, 125320 (2005). 

[21] Th. Martin, R. Landauer, Phys. Rev. B 45, 1742 (1992). 



